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Abstract. We prove that every Markov solution to the three dimen- 
sional Navier-Stokes equation with periodic boundary conditions driven 
by additive Gaussian noise is uniquely ergodic. The convergence to the 
(unique) invariant measure is exponentially fast. 

Moreover, we give a well-posedness criterion for the equations in 
terms of invariant measures. We also analyse the energy balance and 
identify the term which ensures equality in the balance. 



1. Introduction 

The Navier-Stokes equations on the torus with periodic boundary con- 
ditions forced by additive Gaussian noise are a reasonable model for the 
analysis of homogeneous isotropic turbulence for an incompressible New- 
tonian fluid. 

^ (it-vAu + (u-V)u + Vp = 77, 

1 divw = 0. 

The equations share with their deterministic counterpart the well-known 
problems of well-posedness. It is reasonable, and possibly useful, to focus 
on special classes of solutions, having additional properties. 

This paper completes the analysis developed in [11], [12] and [13] (see 
also [ ]). In these papers it was proved that it is possible to show the ex- 
istence of a Markov process which solves the equations. Moreover, under 
some regularity and non-degeneracy assumptions on the co variance of the 
driving noise, it has been shown that the associated Markov transition ker- 
nel is continuous in a space W with a stronger topology (than the topology 
of energy, namely L 2 ) for initial conditions in W . 

In this paper we show that, under suitable regularity assumptions on the 
covariance, every Markov solution admits an invariant measure. Moreover, 
if the noise is non-degenerate, the invariant measure is unique and the con- 
vergence to the (unique) invariant measure is exponentially fast. 
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We stress that similar results have been already obtained by Da Prato & 
Debussche [2], Debussche and Odasso [5] and Odasso [19], for solutions 
obtained as limits of spectral Galerkin approximations to (1.1), and con- 
structed via the Kolmogorov equation associated to the diffusion. The main 
improvement of our results is that such conclusions are generically valid 
for all Markov solutions and not restricted to solutions limit to the Galerkin 
approximations (this would not make any difference whenever the problem 
is well-posed, though) and is general enough to be applied to different prob- 
lems (see for instance [1]). Our analysis is essentially based on the energy 
balance (see Definition 2.4 and Remark 2.5), and in turn shows that such 
balance is the main and crucial ingredient. 

It is worth noticing that the uniquely ergodic results hold for any Markov 
solution, hence different Markov solutions have their own (unique) invariant 
measure. Well-posedness of (1.1) would ensure that the invariant measure 
is unique. We prove that the latter condition is also sufficient, as if only one 
invariant measure is shared among all Markov solutions, then the problem 
is well-posed. 

Finally, we analyse the energy balance for both the process solution to 
the equations and the invariant measure. Due to the lack of regularity of 
trajectories, the energy balance is indeed an inequality. We identify the 
missing term and, under the invariant measure, we relate it to the energy flux 
through wave-numbers. According to both the physical and mathematical 
understanding of the equations, this term should be zero. 

A non-zero compensating term from one side would invalidate the equa- 
tions as a model for phenomenological theories of turbulence, and from the 
other side would show that blow-up is typically true. We stress that neither 
the former nor the latter statements are proved here. 

1.1. Details on results. In the rest of the paper we consider the following 
abstract version of the stochastic Navier-Stokes equations (1.1) above, 

(1.2) du + vAu + B(u,u) = Qj-dW, 

where A is the Laplace operator on the three-dimensional torus T 3 with 
periodic boundary conditions and B is the projection onto the space of 
divergence-free vector fields with finite energy of the Navier-Stokes non- 
linearity (see Section 2.1 for more details). Moreover, W is a cylindrical 
Wiener process on H and a is its covariance operator. We assume that a is 
a symmetric positive operator. We shall need additional assumptions on the 
covariance, as the results contained in the paper holds under slightly dif- 
ferent conditions. Here we gather the different additional assumptions we 
shall use. 

Assumption 1.1. The following assumptions will be used (one at the time) 
throughout the paper. 

[Al] a has finite trace on H. 

3 1 

[A2] there is Oq > such that A* +ai) Q2 is a bounded operator on H. 
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[A3] there is OCq > g such that A* +CXi) Q^ is a bounded operator on H. 

[A4] there is OCq > g swc/z ?/za? A? +a °Q2 is an invertible bounded operator 
on H, with bounded inverse. 

Notice that each of the above conditions implies the following one. We 
shall make clear at every stage of the paper which assumption is used. 

The first main result of the paper concerns the long time behaviour of so- 
lutions to equations (1.2). We show that every Markov solution is uniquely 
ergodic and strongly mixing (Theorem 3.1 and Corollary 3.2). Moreover, 
under an additional technical condition (see Remark 2.5) we prove that the 
convergence to the (unique) invariant measure is exponentially fast (Theo- 
rem 3.3). 

We stress that uniqueness of invariant measure is relative to the Markov 
solution it arises from. As we do not know if the martingale problem associ- 
ated to equations ( 1 .2) is well-posed, in principle there are plenty of Markov 
solutions, and so plenty of invariant measures. In Section 4 we study a few 
properties of the set of invariant measures. In particular, we show the con- 
verse of the above statement, that is if there is only one common invariant 
measure for all Markov solutions, then the martingale problem is well posed 
(Theorem 4.6). 

We also give some remarks on symmetries for the invariant measures 
(such as translations-invariance). Finally, we analyse the energy inequality 
(given as [M3] and [M4] in Definition 2.4, see also Remark 2.5). In particular, 
we identify the missing term in the inequality which, once added, provides 
the equality. For an invariant measure jl, we show that 

1 , 

veGu) + i(jU) = 2° ' 

where \o 2 is the rate of energy injected by the external force, e(ju) = 
E^[|Vx| 2 ] is the mean rate of energy dissipation and is the mean rate 
of inertia! energy dissipation. We show also that is given in terms of 
the energy flux through wave-numbers (see Frisch [15]) as 

l(H) = lim E M £ (x m ■xt){m-xi) 



l+m=k 

\k\-<K, 
\mL>K 
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2. NOTATIONS AND PREVIOUS RESULTS 

2.1. Notations. Let T 3 = [0,2n] 3 and let be the space of infinitely 
differentiable vector fields <p : R 3 — > R 3 that are divergence-free, periodic 
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and 



/ <p(x)dx = 0. 



L 3 

We denote by H the closure of in the norm of L 2 (T 3 ,R 3 ), and similarly 
by V the closure in the norm of H l (T 3 , R 3 ) . Let D(A) be the set of aftueH 
such that AuEH and define the Stokes operator A : D(A) ^HasAu = —Am. 
By properly identifying dual spaces, we have that D(A) C V C H = H' C 
V'cD(A)'. 

The bi-linear operator B : V x V — > V' is defined as 

(5(w,v),w)= V / Wi(x)uj(x) V * X dx 
ifciJTs dxj 

(see Temam [ ] for a more detailed account of all these notations). 

Since A is a linear positive and self-adjoint operator with compact in- 
verse, we can define powers of A. We define two hierarchies of spaces re- 
lated to the problem, using powers of A. The first is given by the following 
spaces of mild regularity (since they are larger than the space V), 

(2.1) v £ =D(A? +£ ), ee(0,~], 

while the second is given by the following spaces of strong regularity, 

(2.2) W a =D(A 6 ^), ae(0,oo), 
where is defined as 



(2.3) 0(a) 



Notice that for every £q and (Xq as above, 




ae (0,i), 



c n> a<) CW = V = ViC V £Q . 

4 

In the proof of most of the results of the paper we shall use repeatedly the 
following inequalities. 

Lemma 2.1 (Temam [22, Lemma 2.1, Part I]). If u G D(A ai ), v G D(A ai ) 
and w G D(A a3 ), then there is a constant cq = co(ofi, Ofe, OC3) such that 

(B(u,v),w) H < c Q \A ai u\ ■ \A a2+ ?v\ ■ |A%v|, 

where a, > and (X\ + «2 + CC3 > 4 if a, 7^ | for all i = 1,2,3, and <X\ + 
OL2 + OC3 > I otherwise. 

Lemma 2.2 ([ , Lemma D.2]). Let a > and u, v G D(A e ^). If a ^ \, 
there is a constant Cq = Cq(cc) such that 

\A a -h(u,v)\ H <C Q \A e ^u\H\A e ^v\H, 
where 6 is the map defined in (2.3). If (X = \, then B maps D(A*) x D(A*) 
continuously to D(A~z~ £ ), for every £ > 0. 
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2.2. Markov solutions to the Navier-Stokes equations. In this section 
we recall a few definitions and result from papers [ ] and [13], with some 
additional remarks. 

2.2.1. Almost sure super-martingales. We say that a process 6 = (9 t ) t >o 
on a probability space (£2,P, J£~), adapted to a filtration (J^) f >o is an a. s. 
super-martingale if it is P-integrable and there is a set T C (0, °°) of null 
Lebesgue measure (that we call the set of exceptional times of 6) such that 

(2.4) E[0 t \^ s ]<0 s , 

for all s £ T and all t > s. 

Lemma 2.3. If 6 = (0 t )t>o is an a. s. super-martingale, then for every s>0 
and every <p e Q°(R) with <p > and Supp<p C [s,°°), 



(2.5) 



E 



(p'(r)0 r dr 



>0. 



Proof. Fix s > and consider a positive smooth map <p with compact sup- 
port in [5,00). By a change of variable, using the a. s. super-martingale 
property, 



E 



((p(r)-(p(r-e))0 r dr 



1 



E 



<p(r)(0 r -e r+£ )rfr 



and in the limit as £ J, we get (2.5). 



>0, 
□ 



It is easy to see that the converse is true (that is, if (2.5) holds, then the 
process is an a. s. super-martingale) under the assumption that the a-fields 
{^t '■ t > 0} are countably generated and is lower semi-continuous (see 
[14]). 

2.2.2. Weak martingale solutions. Let £1 = C([0,°°);D(A)'), let SB be the 
Borel a-field on Q. and let t, : Q. — > 0(A)' be the canonical process on Q. 
(that is, &(fi>) = A filtration can be defined on ^ as ^ = ct(£ s : < 

s<t). 

Definition 2.4. Given \1q e Pr(//), a probability P on (Q.,33) is a solution 
starting at /Jo to the martingale problem associated to the Navier-Stokes 
equations (1.2) if 

[mi] P[Lr oc ([0,oo) ; //)nL 1 2 oc ([0,oo) ; y)] = 1; 

[M2] for each <p e the process Mf , defined P-a. s. on (Q, ^) as 

M? = (&-Zq,<p) h + v [ (% s ,A(p) H ds- [ (B(£ s ,(p),% s ) H ds 

Jo Jo 

is square integrable and (M t ?,3§ t ,P) is a continuous martingale with 

quadratic variation [Af] f = ? | Q, 2 <p |^; 

[M3] the process defined P-a. s. on (Q,^) as 



S 



\&\ 2 H + 2v f \% s \y-ds-tTr[cL) 
Jo 



is P-integrable and 



t ,P) is an a. s. super-martingale; 
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[M4] for each n>2, the process Sf, defined P-a. s. on (CI, 38) as 

£» = \&$+2nv f \^- 2 \md S -n{2n-\)Tv[d] f\^~ 2 ds 
Jo Jo 

is P-integrable and (S^ 1 ,38 t ,P) is an a. s. super-martingale; 

[M5] jUo is the marginal of P at time t = 0. 

Remark 2.5 (enhanced martingale solutions). A slightly different approach 
has been followed in [1] to show existence of Markov solution for a different 
model (an equation for surface growth driven by space-time white noise), 
as the energy balance has been given in terms of an almost sure property. In 
the Navier-Stokes setting of this paper the property reads (some equivalent 
statements are possible as in [1]) 

[M3-as] there is a set Tp x C (0,°°) of null Lebesgue measure such that for all 
s £ Tp x and all t>s, 

W(v,z)<^(v,z)] = l, 

where Sf is defined as 

1 f 1 f 1 

&t(v,z) = -\v t \ 2 H + V \v r \ldr+ (v r ,B(v r + z r ,z r ))Hdr, 

2 Jo Jo 

z is the solution to the Stokes problem (A. 2) and v = £ — z. It is possible to 
show that, as in [ ], there exist Markov solutions which additionally satisfy 
[M3-as]. We shall assume this statement (see [14] for more details). 

2.3. Previous results. In the next theorems we summarise some results on 
existence and regularity of Markov solutions to the Navier-Stokes equations 
(1.2). First we show that there is a Markov solution to the Navier-Stokes 
equations (1.2). 

Theorem 2.6 ([12, Theorem 4.1]). Under condition [Al] of Assumption 1.1, 
there exists a family (P x ) x ^h of weak martingale solutions (as defined above 
in Definition 2.4), with P x starting at the measure concentrated on x, for 
eachx G H, and the almost sure Markov property holds. More precisely, for 
every xEH there is asetTc (0, °°) of null Lebesgue measure such that for 
all s^T,allt>s and all bounded measurable (j) : H — > R, 

E^[0(^)|^]=E^[0(^)]. 

The map x i— > P x is in principle, from the above result, only measurable. 
The regularity of dependence from initial condition can be significantly im- 
proved under stronger assumptions on the noise, as shown by the theorem 
below. 

If (P x )xeH is a Markov solution, the transition semigroup 1 associated to 
the solution is defined as 

(2.6) 2>,<p(*)=E p *[<p(&)], xeH,t>0, 

^Notice that, due to the Markov property holding only almost surely, the family of 
operators (2>r) r >o is n °t a semigroup, as the semigroup property holds for almost every 
time. 
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for every bounded measurable (p : H ^ R. 

Theorem 2.7 ([12, Theorem 5.11]). Under condition [A4] of Assumption 
1.1, the transition semigroup (fPt)t>o associated to every Markov solution 
(Px)xeH is strong Feller in the topology of W^. More precisely, tp t (j> e 
Cbi^ao) far every t > and every bounded measurable : H — > R. 

The regularity result can be given more explicitly in terms of quasi- 
Lipschitz regularity (that is, Lipschitz up to a logarithmic correction) as in 
[13], albeit the estimate given there holds true only for (Xq = | (an extension 
to all values of (Xq > \ can be found in [14]). 

3. Existence and uniqueness of the invariant measure 

In this section we prove existence of invariant measures by means of the 
classical Krylov-Bogoliubov method. Let (P x ) x eH be a Markov solution 
and denote by (v t )t>o its transition semigroup (see (2.6)). Let xq<eH and 



where is the Dirac measure concentrated on xq. It is known (see for 
example Da Prato & Zabczyk [ ]) that any limit point of the family of prob- 
ability measures (Ht)t>o is an invariant measure for (2v) f >o, provided that 
the family is tight in the topology where the transition semigroup is Feller. 

Theorem 3.1. Assume [A2] of Assumption 1.1. Let (P x ) x eH be any Markov 
solution to the Navier-Stokes equations (see Theorem 2.6) and let (p t )t>0 
be the associated transition semigroup. Then the family of probability mea- 
sures (Ht)t>l is tight in Wo®- 

The above theorem, together with the strong Feller property ensured by 
Theorem 2.7 and Doob's theorem (see Da Prato & Zabczyk [ ]), immedi- 
ately imply the following result. 

Corollary 3.2. Under [A4] of Assumption 1.1, every Markov selection to 
the Navier-Stokes equations has a unique invariant measure /A*, which is 
strongly mixing. Moreover, there are 8 > and y > (depending only on 
(Xq) such that 

The convergence of transition probabilities to the unique invariant mea- 
sure can be further improved if, under the same assumptions of above re- 
sults, we deal with the enhanced martingale solutions introduced in Remark 
2.5. This is a technical requirement that makes the proof of Theorem 3.3 
below simple and, above all, feasible. 

Theorem 3.3. Assume [A4] of Assumption 1.1 and consider an arbitrary 
Markov solution (P x ) x ^h made of enhanced martingale solutions (see Re- 
mark 2.5). Let /l* be its unique invariant measure. Then there are constants 



(3.1) 
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C exp > and a > ( independent of the Markov solution and depending only 
on the data of problem) such that 

\\t*8 Xq - TV < C exp (\ + \xo\ 2 H )e~ at , 

for all t > and xq G H, where \\ ■ \\ TV is the total variation distance on 
measures. 

Remark 3.4. The proof of Theorem 3.3 above actually shows a slightly 
stronger convergence, namely 



sup |av0(x o )- / ^>(x)n(dx)\ < C exp {\ + \x\ 2 H )e 
W\v<i J 



at 



for every x E H and t > 0, with same constants C exp and a, where the norm 
|| ■ || v is defined on Borel measurable maps § : H — > R as 

ii , I, 10 (*)| 

xeH 1 + \X\fj 

(see Goldys & Maslowski [ ] for details). 

From Theorem 13 of [13] and again from Theorem 4.2. 1 of Da Prato & 
Zabczyk [ ] we also deduce the following result. 

Corollary 3.5. Under the assumptions of previous corollary, let \l\ and /l2 

be the invariant measures associated to two different Markov selections. 
Then the two measures are mutually equivalent. 

The rest of the section is devoted to the proof of Theorems 3.1 and 3.3. 

3.1. The proof of Theorem 3.1. We fix a Markov solution (P x ) x ^h- Prior 
to the proof of the theorem, we show two lemmas on momenta of the solu- 
tion. The second lemma is the crucial one. 

Lemma 3.6. Assume [Ai] of Assumption 1.1. Then for every x G H and 
t>0, 

^[I^I^^W^ + ^ll-e- 2 -). 

Proof. The result easily follows from the super-martingale property [M3], 
Poincare inequality and Gronwall's lemma (see for example [20] for de- 
tails). □ 

Lemma 3.7. Assume [A2] of Assumption 1.1. Then there are C > 0, 8 > 
and 7 > depending only on £q, (Xq, v and o 2 (but not on the Markov 
solution) such that for xq G H and t > 1, 



f\A 5 ^ ds]<C(l + \x \ 2 H ) 
t Jo "o 



A slight modification of the argument in the proof below provides an 
inequality similar to that of the lemma also for t < 1 . 
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Proof. Let £o G (0, |] with £o < 2«o- We first prove the statement of the 
lemma for xq G •U £o . 

Consider values 8 = 5(eo, (Xo), y = y(£o, olq) provided by Theorem A. 2. 
For every fixed value M > we choose R > 1 + 2|jco|#> whose value will be 
given explicitly later, and we denote by £r the small time where the blow-up 
estimate (A. 5) of Theorem A.l holds true. 

Fix t > 1 and £ < e R , and let n £ G N be the largest integer such that 
£n £ < t. By the Markov property, 



j ie rke+e 



(3.2) 



< ? 

= -YE' 
t h n 



P*[\A*S,\\ >M]ds 



rke+e 



P^[\A 5 ^_ k£ \l % >M]ds 



1 " E r i 



" k=0 

where /l f is the measure defined in (3.1). Now, by Theorem A.l, for every 
x G l^en such that Ixl?, < R, 

and so, by using (A. 5) and Chebychev inequality, 

f&[|A a i;i 2 ^>M]< 

< (ifllA'i'li, > M] +P {l ,W-'> < , Ri + , 

n) 



<h\F fl ^r\ + — E ^ [|A 5 £ v ry ]+cie" C2 ^. 
— listed >«} ^/y u ^^oq 1 1 

We use the above inequality in (3.2) and we apply Theorem A. 2 and the 
previous lemma, 

Vt[\A S x\ 2 >M]< 



1 ,l£ r 1 P ( £ o- R ) r£ . « 2 n 



1 p(«b,«) re 



< jl^O^i 2 ^ >R]+j^(l + e+^[\^\%]) + c 1 ec^ 
F « £ r r 2 



< 



< 
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Since all computations above are true for all £ < e^, if we integrate for 
£ G (J;£fl,£fl), we get 



< 



£r 



f k=0 



Mr 



We use the energy inequality and the previous lemma to estimate the only 
complicated term in the inequality above, 

' J t k=l * i— 1 •/-*- 



k=\ ' 

Fa " R 

< ^ y e p *o 



k=l 
nR j 



Br 



ke R 



\%r\v dr 



<^log^-(l + h||), 

rt £/? 

where we have set = Ha and n e < for all £ G [|cr, £#]. Since by (A.5) 
the dependence of £^ on R is like i? a , for some exponent a depending on 
£o, we may choose R in such a way that for every t > 1, 

Alr[|A**|2 > M] < J_i ogM . 



1 



M b 



for a suitable > 0. In conclusion, the statement of the lemma is proved for 
initial conditions xq G n/^. 

If xq G H, since for every s>0we know that E, s G i4 , P^-a. s., then by 
the Markov property, 



E^o [ / ^ft rf r ] = E^o [E^. [ / \A S C\ 2 " dr}] 



t<5e/|2y 



l «0 



<a(i + |x |ff), 

where we have used the previous lemma and this same lemma for initial 
conditions in ^ £o . Finally, as s [ 0, the conclusion follows by the monotone 
convergence theorem. □ 

Proof of Theorem 3.1. Choose an arbitrary point xq G H and consider the 
sequence of measures (fJ.t)t>l defined by formula (3.1). Since 

/iAv w >(&)=}/ o 'iA^e ao *, 
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where the constants 8 and y are those provided by the previous lemma, 
it follows by that same lemma that the sequence of measures is tight in 

WW □ 

3.2. The proof of Theorem 3.3. As stated in the statement of the theorem, 
in this section we work with the enhanced martingale solutions defined 
in Remark 2.5. It means that the energy balance [M3-as] is available for 
proofs. Prior to the proof, we give a few auxiliary results, summarised in 
the following lemmas. In the first one we show that any solution enters in a 
ball of small energy with positive probability. 

Lemma 3.8 (entrance time in a ball of small energy). Assume [A3]. Given 
R > and 8 > 0, there exists T\ = T\(8,R) such that 

inf P(T h x,{y:\y\ 2 H <8})>0. 

\ X \ 2 H <R 

Proof. Consider a value k\ = k\{8), to be chosen later, and let A = {ft) : 
su Pf6[0,7i] \ z t\H — ^i}- We know that for every x, the value P X [A] > is 
constant with respect to x. Since \^\h < \zt\H + \ v t\H, we shall estimate v. 
For all ft) G A such that the inequality in [M3-as] (at page 6) holds, we have 

\ v i\h - \v s \h + 2v f \v r \ydr<cf \z r \v\^ r \v\A^v r \ H dr 

J s J s 

ft I 1 I 3 

(3-3) - C J. (\ Zr \ V \ Vr \v\ V r\v + \ Z r\v\ V r\H\Vr\v)dr 

- y J s \ v r\vdr + k 2 J^ \v r \ 2 H dr + k 3 (t-s), 

A - - 

where we have set k 2 = c(k\ + k\ ) and kj = ck\ . By the Poincare inequality 
(the first eigenvalue of the Laplace operator on the torus T 3 is 1), 

\v t \ 2 H + ( v - k 2) j \v r \ 2 H dr< \v s \ 2 H + k 3 (t-s), 
and Gronwall's lemma ensures that 



\l < \x\W {v ~ h)t + -^V(l - z- {v - kl)t ) < R^-W + • 



V - k 2 ~~ V - k 2 

If we choose k\ and T\ in such a way that 

ki<jS, k 2 <v, -^r<k Re^-W<\8, 
4 v-k 2 8 8 

we finally obtain that, if |jc|| < R, then Px[{\&i \h < 5 ) nA ] = p x[ A ]- n 

The second lemma shows that with positive probability the dynamics en- 
ters into a (sufficiently large) ball of space V. 
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Lemma 3.9 (entrance time in a ball of finite dissipation). Assume [A3] from 
Assumption 1.1. Then there exists 8 > small enough such that there are 
T 2 = T 2 (S) > andR 2 = R 2 (S) > with 

inf P(T 2 ,x, {y:\y\l <R 2 })>0. 
\x\ 2 H <5 

Proof. Set T 2 = 1 and let A = {supr 01 i |A§z|^ < k{\, with k\ to be chosen 
later, together with 8. 

For all ft) G A for which the inequality in [M3-as] (at page 6) holds, we can 
proceed as in the proof of Lemma 3.8 to get \v t \jj < 8 + ^d^, with k 2 and 
£3 defined similarly. Using (3.3), we get 

/ \v s \yds<- 8+k 3 +k 2 (8 + — ) 
Jo v L v v - k 2 ' 

where k\ is small enough so that k 2 < v. 

Next, we notice that the set {r 6 [0, 1] : |v r |y < 2^4} is non-empty (its 
Lebesgue measure is larger than one half). So for each ro in such a set, 
|vr lv — ^^4- Since the energy inequality [M3-as] holds, for a short time after 
ro, v coincides with the unique regular solution. We shall choose k\ and 8 
small enough so that the short time goes well beyond 1 . 

Indeed, using (2.1) (as in (A. 6) with £0 = |), we get for suitable universal 
constants c\ and c 2 , 

^|v|^ + 2v|Av|^<v|Av|^ + c 1 (|v|^ + |Ai z | 4 ^)<v|Av|^ + c 1 (|v|^ + c2fc 1 ), 

- - 
and so, if <p(r) = \v r \y + k%, we have <p(ro) < 2k$ + k\ and <p < c\(p 3 . Now, 

if we choose k\ and 8 small enough so that 

Aci{lk A + k\) 2 < 1 

the solution to the differential inequality of <p is finite at least up to time 
1 -fro. In particular, <p(l) < (2ci)~2 and so by easy computations, 

2 



\& 2 \ v = l&lv < (|vi|v + |zi|v) 2 < 2h + 



We choose now the last term on the right-hand side of the above formula as 
R 2 . In conclusion, P(T 2 ,x, {\y\ v < R2) > Px[A] and again the value of P X [A] 
is independent of x. □ 

In the last auxiliary lemma we show that the dynamics enters in a compact 
subset of tVciq. This is crucial since the strong Feller property holds in the 
topology of "jVo,, (Theorem 2.7). 

Lemma 3.10 (entrance time in a ball of high regularity). Assume [A3] from 
Assumption 1.1. Then there is /3 > (depending only on (Xq) such that for 
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every R2 > there are a time T 3 = T 3 (R2) > and a constant C = C(/?2) > 
and 

inf P(T 3 ,x,{y:\APy\l < C}) > 0. 

Proof. Given i?2 > 0, we choose /3 = 0" — 9((Xq), T 3 and C as given in 
Lemma A. 3. Notice that the set K = {y : |A^y|^, < C} is a compact subset 

of tVoq. 

If T = t (j ' 3R2) is the time up to which all solutions starting at x coincide 
with the unique solution to problem (A. 1), then 

P(T 3 ,x,K)=P x [\A^T 3 \l ao <C) 

>P,[\A^ T3 \l ao <c, r>r 3 ] 

>P^ ] [\AP^\l, aa <C,T>T 3 ]. 

Now, the conclusion follows from Lemma A. 3. □ 

Proof of Theorem 3.3. Let (P x ) x eH be a Markov solution and consider the 
corresponding transition kernel (P(t,x, -))t>0jceH- We choose the value of 
Co given in Lemma A. 3 and we consider the value 6" > 0(ab) provided by 
the same lemma. 

The exponential convergence follows from an abstract result of Goldys 
& Maslowski [16, Theorem 3.1] (which, in turns, is based on results from 
the book by Meyn & Tweedie [17]). More precisely, we need to verify the 
following four conditions, 

1. the measures (P(t, x, -))t>o^efl 316 equivalent, 

2. x — > P(t,x, r) is continuous in for all t > and Borel sets T C H, 

3. For each R > 1 there are Tq > and a compact subset £ C such 
that 

inf P(T ,x,K) >0, 

4. there are k, b, c > such that for all t > 0, 

E*[|&&]<*M|e-* + c. 

The first property follows from Theorem 13 in [13] (there equivalence is 
stated only for x E tVoq, but it easy to see by the Markov property that it 
holds for x E H, as •W ao is a set of full measure for each P(t,x, •)). The 
second property follows from the strong Feller property, while the fourth 
property follows from Lemma 3.6. 

We only need to prove the third property. We fix R > 1 and we wish to 
prove that there are Tq = Tq(R) and K = K(R) such that 

(3.4) inf P(T ,x, {y : \A°"y\ 2 H < K}) > 0. 

\x\ 2 H <R 

We choose the value 8 provided by Lemma 3.9 together with the time T2 
and value R2. Corresponding to the values R and 8, Lemma 3.8 gives a time 
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T\. Moreover, corresponding to R 2 , Lemma 3.10 provides the time T 3 and 
value C. 

We set Tq = 7i + T 2 + ?3, then if \x\g < R, using three times the Markov 
property, 



P(T ,x,K) > 

> inf P(T 3 ,z,K) inf P{T 2 ,y, {\z\ 2 v < R 2 }) inf P(T u x,{\y\ 2 H < 8}) 



and the right-hand side is positive (and bounded from below independently 
of x) due to Lemma 3.8, 3.9 and 3.10. 

Finally, the constants C exp and a in the statement of the theorem are inde- 
pendent of the Markov solution since all computations either depend on the 
data (the viscosity v, the strength of the noise a 2 , etc., such as in Lemma 
3.6) or are made on the regularised problem analysed in the appendix. □ 



In the previous section we have shown that, under suitable assumptions 
on the driving noise, every Markov solution has a unique invariant mea- 
sure. As in principle there can be several different Markov solutions, so are 
invariant measures. 

In the first part of the section we show that well-posedness of the martin- 
gale problem associated to (1.2) is equivalent to the statement that there is 
only one invariant measure, regardless of the multiplicity of solutions. 

In the second part we give some remarks on symmetries of invariant mea- 
sures, while in the third part we analyse the energy balance. 

4.1. A connection between uniqueness of invariant measures and well- 
posedness of the martingale problem. 

4.1.1. Stationary solutions. Consider the (unique) invariant measure asso- 
ciated to a Markov solution (P x )xeH, as provided by Corollary 3.2, and de- 
fine the following probability measure 



Lemma 4.1. The probability measure P+ defined above is invariant ( in the 
following, stationary) with respect to the time shifts r\ t :Q.^ Q. defined as 



Proof. It is sufficient to prove that the finite dimensional marginals of P+ and 
T] S P* are the same. The case of one single time is easy, by invariance of fu. 
We consider only the two-dimensional case (one can proceed by induction 
in the general case). Consider t\ <t 2 , then by the Markov property and 



4. Further analysis of equilibrium states 



(4.1) 




TJf(fl>)(j) = G)(t + s). 
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invariance of jLi*, 

E^[/(^ 2 )] = / E p *[M s+tl ,Ut 2 )}Lh(dx) 



e p x [E p ^n){f(eo( S + h), ^ )]J 
E P *[F(&W] = | E^'[F(^)]au(^) 
E^[/(^,4)]^)=E^[/(4,^)], 

where in the above formula we have setF(j) = W>'[f(y, <^ 2 _ fl )]. □ 

In turns, the lemma above ensures that is the unique probability mea- 
sure on Q. such that 

1. P± is stationary, 

2. is associated 2 to the Markov solution (P x )x<eh- 

Uniqueness follows easily since [l* is the unique invariant measure of the 
Markov solution (P x ) x <eh an d since the law of a Markov process is deter- 
mined by its one-dimensional (with respect to time) marginal distributions 
(as in the proof of Lemma above). We shall see later on that for a special 
class of invariant measures this uniqueness statement can be strengthened 
(see Proposition 4.4). 

In general one can have several stationary solutions (see for example [20] 
for the definition and a different proof of existence) and possibly not all of 
them are associated to a Markov solution. Hence we define the two sets, 

(4.2) / = {jiie Pr(H) : /i is the marginal of a stationary solution} 

(4 3) J? = ill e Pr(H) • ^ 1S tne urnc l ue invariant measure associated 1 
m V ^ ' ' to a Markov solution J 

and, trivially, J m C J . 

Remark 4.2 (Topological properties of and J^ m ). By the same properties 
that ensure existence of solutions (and following similar computations, see 
for example [ ]), it is easy to see that is a compact subset of £1. More- 
over, by Corollary 3.2, <0 m and hence J* e are relatively compact in a much 
stronger topology. 

4.1.2. A short recap on the selection principle. It is necessary to give a 
short account on the procedure which proves the existence of Markov se- 
lection (namely, the proof of Theorem 2.6). We refer to [12] for all details. 

Given x E H, let ^(x) C Pr(Q) be the set of all weak martingale solutions 
(according to Definition 2.4) to equation (1.2), starting at*. 

In the proof of Theorem 2.6 (see [ ]) the sets ^(x) are shrunken to one 
single element in the following way. Fix a family f n ) n > l which is dense 



We say that a probability measure P on Q. is associated to a Markov solution {P x )xeH 
if for every t > 0, P\% = Pm{t) for P-a. e. co S i2, where (P\% )<oe£2 is a regular conditional 
probability distribution of P given SB t . 
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in [0,°o) x Cb(D(A)') and consider the functionals J n = Jx n j n , where 3% * is 
given by 



r r°° 

J l j(P)=E p [j Q c- Xt f(^ t )dt . 



for arbitrary A > and / : D(A)' — > R upper semi-continuous. Next, set 



^(x), % l (x) = {PeV n -i{x):Jn(P)= sup J n (Q)}. 



All these sets are compact and their intersection is a single element (the 
selection associated to this maximised sequence), DneN^W = {Px}- 

Given now a probability measure /i on D(A)', one can define the set < «f (ji) 
as the set of all probability measures P on £1 such that 

1. the marginal at time of P is /i; 

2. there is a a map jc 1 — > Q x : H — > Pr(Q) such that P = J Q x \i{dx) and 
<2 X G ^(x) for all x (in different words, the conditional distribution 
of P at time is made of elements from sets ( < t£(x)) X £h). 

We can now give the following extension to the selection principle. 

Proposition 4.3. Let (P x ) x eH be the Markov selection associated to the se- 
quence (A„, f n )n>h Then the probability P^ = J H P x fi(dx) is the unique 
maximiser associated to the sequence {Kiifn)ri>i- More precisely, 



Proof. Since each Q G ^(/i) is given by Q = j Q x /i(dx), for some x t— > Q x , 
by linearity of the map 7i it easily follows that P^ G ^i(ju). Moreover, 
each Q G ^1 (/i) has a similar structure: Q = J Q x /J, (dx) and <2 X G ^ (x) for 
/i-a. e. x e H. In fact, /i(<2x) < Ji[P x ), M _a - s -, an d -A(Q) = -M^u)> an d 
so /i(<2jc) = Ji(P.x), for /i-a. e. x. By induction, G ^(jii) and for each 
2 = jQxH(dx) G Q x G <*f„(x), for /i-a. e.xeH. 

In conclusion, G 'C(m) = fl^M) and for each Q = / Q x [i(dx) G 
(/i), <2x G ^(x), for /i-a. e. x £ H. But we know that each ^L(x) has 
exactly one element, P x , so that in conclusion the only element of ^(/i) is 



4.1.3. Connection with well-posedness. Now, if we are given a sequence 
(hnfn)neti as above, the selection principle provides a Markov solution 
(Px)xeH- Corollary 3.2 ensures that this Markov solution has a unique in- 
variant measure /i*. Moreover, from the proposition above, the stationary 
solution P* = J i> x /i* (dx) is the unique sequential maximiser of the sequence 



sup Ji(P), 




sup J„(P), 
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(•Ai)neN on ^Qtl*). This justifies, in analogy with the definition of (4.2) and 
(4.3), the definition of the following set, 

{H is the invariant measure associated to a Markovs 
\l : solution obtained by the maximisation procedure for > , 
some sequence (A„,/„)„ eN J 

and obviously C J^ ra C J' . 

Proposition 4.4. If jl* G y e , then the stationary solution associated to 
jU* is the unique stationary measure in ^(/i*). 

Proof. Since /l* G by definition there is a sequence (A„,/ m )„ g n dense 
in [0,°°) x Cb(D(A)') such that P* maximises functionals J n = Jx n j„ (one 

after the other, as explained in Proposition 4.3). Now, if P G ^(ju*) is a 
stationary solution, then 

J n {P) =E p [j ° Q -^ t f n ^ t )dt\ = (| f n (x)n*(dx)) [z~ Kt dt 

= J- J fn(x)H*(dx), 

and so J n (P) = J n (P*) for all n. By Proposition 4.3, it follows that P = 
P*. □ 

If we consider now Markov solutions as those obtained for the Navier- 
Stokes equations, namely each of them is strong Feller and irreducible on 
t^oq, the previous result gives immediately a criterion for well-posedness. 
In few words, uniqueness of the invariant measures among Markov solu- 
tions is equivalent to well-posedness of the martingale problem. 

Corollary 4.5. Assume that every Markov selection is W^-strong Feller 
and fully supported on W^. If(Px)xeH and (P x ) x€ h are two Markov selec- 
tions, with (P x ) x eH coming from a maximisation procedure, and they have 
the same invariant measure, then they coincide on "Woq. 

Proof. Let P+ and P{ be the stationary solutions associated to the two selec- 
tions. If the two selections have the same invariant measure, it follows from 
the previous theorem that they have the same stationary solution, that is 
P+ = P+. It follows from this that their conditional probability distributions 
at time coincide, 

P X = P' X fx* — almost surely. 

By W^-strong Feller regularity and irreducibility they coincide on every 

*G^V □ 

We summarise the result in the following theorem. It follows easily from 
the previous corollary and from the fact that well-posedness of the martin- 
gale problem is equivalent to uniqueness of Markov selections (see Theo- 
rem 12.2.4 of Stroock & Varadhan [21]). 
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Theorem 4.6. Under [A4] of Assumption 1.1, assume that the set J? e defined 
in (4.4) contains only one invariant measure. Then the martingale problem 
associated to the Navier-Stokes equations (1.2) is well-posed on 'W^ {and 
hence on "V^for all Co > 0). 

Proof. We only have to prove that, given two Markov solutions (P x ) x eH and 
{P' x )x<eh, for every x E we have P x = P' x . This statement holds for x E fVa 
by the previous corollary. Fix Co > and x E n/^ . Choose R > > \x\ r, , then, 
for every bounded continuous 0, by the Markov property, 

¥ t (j)(x) = E p *[v t _ s (i)(^_ s )l {r{£o , R]>§} } +E p *[¥ t _ s (j)(% t _ § )t {T{£o , R)>§} } 

where v is the transition semigroup of (P x )xeH- The first term on the right- 
hand side is independent of the selection, by the weak-strong uniqueness of 
Theorem A. 1 , hence 

$ t Q{x)-<p' t $(x) = 

= E p *[¥ t _ s <i>(Z t _ s )t { ^ R)>S} ] +E%'_ g 0(4_ 5 )l {T(8b , fi)>g} ] 

and, by the blow-up estimate of Theorem A.l, as 8 — > 0, we get i>/0(;e) = 
(P t (j)(x) for all <j) and all t. □ 

4.2. Translations-invariance and other symmetries. In the analysis of 
homogeneous isotropic turbulence, for which equations (1.1) can be con- 
sidered a model, it is interesting to consider equilibrium states invariant 
with respect to several symmetries (see for example Frisch [15]). 

Here we are interested in solutions which are translations-invariant (in 
the physical space). For every a E R 3 , define on the map m a : — > 
as 

m a ((p)(x) = cp(a + x), xER 3 
for any (p E 3°° . The map obviously extends to H and D(A a ) for each 
a. By composition, it extends to continuous functions on H (or D[A a ) for 
every a) and, by duality, to probability measures on H. It also extends to Q. 
as 

m a (C0)(t)=m a ((O(t)), t>0,Q)EQ., 
and, by duality, to probability measures on Q.. 

A function (or a measure) is translations -invariant if it is invariant un- 
der the action of (m a ) aieR 3. The Navier-Stokes equations are translations- 
invariant, so equation (1.2) is translations-invariant only if such is the noise. 
The driving noise is translations-invariant if and only if the covariance Q, 
commutes with all m a . It is easy to verify that this is equivalent to have 
homogeneous noise which, in turns, is equivalent to have that Q, and A 
commute. So, easy examples of homogeneous noise compatible with the 

3 

properties of Assumption 1.1 are Q, =A~2~ (X{> for any Oq in the correct 
range. 

Proposition 4.7. Assume that Q, and A commute. Then the following prop- 
erties hold true. 
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1. For every a E R 3 , m a is a one-to-one map on JP, on J? m and on J? e . 

2. There is at least one translations-invariant measure in J? . 

Proof. We first show that if P is the law of a solution to equations (1.2), 
then m a P is also a solution for every a E R 3 . Since for every a E R 3 the 
map m a is an isometry on H, the image of a cylindrical Wiener process on 
H is again a cylindrical Wiener process. The assumption on q, ensures now 
that the noise term is translations-invariant and so it is easy to check that 
all requirements of either Definition 2.4 or of any definition of solutions for 
the stochastic PDE (1.1) available in the literature (see for example Flandoli 
& Gatarek [ ], we refer also to [20] as regarding stationary solutions) are 
verified. 

In particular, if P is stationary, then m a P is again stationary and so m a 
is a one-to-one map on J* . Moreover, since is closed and convex (see 
Remark 4.2), it follows that there exists a translations-invariant measure. 
Indeed, given fl E J?, there is a stationary solution P^ whose marginal is /i. 
Now, the probability measure 

P = t=— r / m„Pn da 
|T 3 | Jt 3 m 

is again a stationary solution and its marginal is translations-invariant, as 
m a +2xk = m a for every k E Z 3 . 

We next prove that m a maps one-to-one. Let /I* E J* m and consider 
a Markov solution (P x )x<eh having fu as one of its invariant measures. Fix 
a E R 3 and set Q x = m a (P m _ a ^). It is easy to verify that (Q x )x<eh is another 
Markov solution, since 

Qx I %, = m a {Pm-a (x) )\%,= m a (^m_ a (ffl) (t) ) P m- a (x) ~ a - s ■ 

Moreover, m a (ju) is an invariant measure of (Q x )xeH- 

Finally, in order to show that m a maps one-to-one, we only need to 
find a maximising sequence for the solution (Q x )xeH defined above. Let 
(^*,/»)neN be a maximising sequence for (P x ) x <eh, then (A„,/„ o m_ a ),, eN 
is a maximising sequence for (Q x )xsh- D 

We stress that in the proposition above existence of a translations invari- 
ant equilibrium measure is granted in y, but we do not know if such a 
measure belongs to y m . 

Notice finally that if problem (1.2) is well-posed, it follows easily that 
the unique invariant measure must be translations-invariant. 

Similar conclusions can be found for other symmetries of the torus, such 
as isotropy (invariance with respect to rotations, see for example [ ] where 
such symmetries are discussed in view of a connections between homoge- 
neous turbulence and equations (1.1)). 
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4.3. The balance of energy. In the framework of Markov solutions ex- 
amined in this paper, the balance of energy corresponds to the a. s. super- 
martingale property [M3] (and, more generally, of [M4]) of Definition 2.4. As 
clarified in [12], the two facts 

1. the balance holds only for almost every time, 

2. the balance is an inequality, rather than an equality, 

correspond to a lack of regularity, in time in the first case and in space in 
the second, of solutions to the equations (1.1). From the point of view of 
the model, such facts translate to a loss of energy in the balance. 

Generally speaking, the problem could be approached by using the Doob- 
Meyer decomposition (which may hold even in this case, where the energy- 
balance process <§ x is not continuous and the filtration (SS t )t>o does not 
satisfy the usual conditions, see Dellacherie & Meyer [6]). We shall follow 
a different approach, due to the lack of regularity of trajectories solutions to 
the equations. We shall see that the bounded variation term in the decom- 
position of S 1 is a distribution valued process. 

Let a > and define the operator L a = exp(— aAz). Given a martingale 
solution P x starting at some x E H, there is a Wiener process (W t ) t >o such 
that the canonical process ^ on Q. solves (1.2). The process L a £, under P x is 
regular enough so that we can use the standard stochastic calculus. Given 
an arbitrary <p e Q°(R), with support in [0,°°), Ito formula gives 



d[<p{t)\L a $ t \ 2 H ] = (p'(t)\L a ^ t \ 2 H + 2(p(t)(L a ^,dL a ^) H 

= <p'(t) \L a $ t || - 2v<p(f ) - 2fl»(0 (L a &,L a B(Z t , $)) H 



P x -Si. s. As a i 0, the operator L a approximates the identity, so that by the 
regularity of ^ under P x , 



+ 29(0 ( L 2a&, d 7 -dW) H + (p(t)a 2 , 



where 



(7% — Tr[o,L2a], and so, by integrating in time, 
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P x -a. s. and in L (Q,), where o 2 = Tr[oJ. In conclusion, the limit 

(4.5) / Jf t (?;)(p(t)dt:=)im f°° (p(t)(L a ^,L a B^ t ^ t ))dt 

J a[0 JO 

exists P x -a. s. and in L 1 (Q.), and defines a distributions-valued random vari- 
able. Moreover, J? depends only on ^ (that is, on P x ) and not on the 
approximation operators (L a ) a> o used. We finally have 

(4.6) 2v J <p(t)\$\ldt + 2 J <p{t)f t {t;)dt = 
(p'{t)\5 )t \ 1 H dt + 2 J (p(t)($ t ,(i?dW)Hdt + o 2 J (p(t)dt, 



The previous computations and Lemma 2.3 provide finally the following 
result. In few words, the next theorem states that the term ^(^) plays the 
role of the increasing process in the Doob-Meyer decomposition of the a. s. 
super-martingale [S t )t>o (defined by property [M3] of Definition 2.4). 

Theorem 4.8. Given a martingale solution P x , there exists a distribution- 
valued random variable ^ defined by (4.5), such that the (distribution- 
valued) process ${ + 2E Px [J J? r {£,) dr\SS t ] is a distribution-valued martin- 
gale, that is for every s>0 and every <p e Q°(R) with Supp<p C [s,°°), 



<pf(t)(4+m p *[ j r {^)dr\m^)dt 



o. 



Moreover, ^ is a positive distribution, in the sense that for every <p e 
C"(R) with <p > and Supp<p C [5,00), 



(4.7) 



E- 



p, 



>0. 



(p{r)f r {t;)dri 
Proof. The first part of the theorem follows easily, since 
J <$>\t)S}dt = J <p?(t)(\$\l + 2v J \$r\vdr-o 2 ?) dt 

= J<p'(t)\$\%dt-2v J ' <p(t)\^dt + o 2 J <p(t)dt, 

and so, using the above computation and formula (4.6), we get the con- 
clusion. The second part is a consequence of the first part (the martingale 
property) and the fact that [S t );>o is an a. s. super-martingale. □ 

Remark 4.9. The Ito formula applied to ^(t)\L a ^\jj provides an analysis 
of the a. s. super-martingale <§ n , defined in property [M4] of Definition 2.4, 
similar to that developed above for S l and ^ (<§ ) . 

Remark 4.10. Duchon & Robert [ ] show that the energy equality holds 
for suitable weak solutions to the deterministic Navier-Stokes equations if 
one takes into account the additional term 9i, a distribution in space and 
time, obtained by means of the limit of space-time regularisations. Their 
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computations in our setting would lead to a random distribution @(%)(t,x) 
in space and time and 

jm)=muT 3 )= I Sim,*)*. 

This is only formal because in principle our solutions are not suitable (see 
[20] for existence of suitable solutions in the stochastic setting). 

Moreover, they relate the quantity Qi to the four-fifth law in turbulence 
theory (see for instance Frisch [ 1 5]). 

4.3 . 1 . The mean rate ofinertial energy dissipation. Consider a Markov so- 
lution (P x )xeH and let fu be its unique invariant measure. Define the mean 
rate of energy dissipation as 

e(^):=E^[|x|2]. 

We know that 2ve( i u*) < o 2 . We can as well consider the expectation with 
respect to the stationary solution of the distribution J? defined in the 
previous section. As /1* is an invariant measure, the distribution <p — > 
E^* [{^/(^), <p)) is invariant with respect to time-shifts. Hence there is a 
constant t(/i*), that we call mean rate ofinertial energy dissipation, such 
that 

^ [Jmz MO dt \ = 1 W f <p (0 dx - 

We notice that, as a consequence of (4.7), 

0<i(M*). 

By taking the expectation in the balance of energy given in (4.6), we 
finally obtain the following energy equality, 

(4.8) ve(^) + t(^) = ^ 2 - 

The quantity can be given as the expectation of (4.5). Notice that in 
this case the expectation in /1* and the limit in (4.5) commute. 

We give a different formulation of t(/l*) in terms of Fourier modes. As 
the definition of Jf (and hence of t) is independent of the approximation 
(as long as the approximating quantities are regular enough, so that all the 
computations are correct), we use a ultraviolet cut-off in the Fourier space. 
For every threshold K, define the projection P^ of H onto low modes as 

P l K x= £ x k e h \ forx= £x,e fc , 
\k\«,<N k£l? 

and the projection onto high modes P^ = / — P^. Applying Ito formula on 
<)!>(£) |P^| 2 as in the previous section, taking the expectation with respect 
to Py^ and then getting the limit as K j 00 yields the following representation 
formula for 

i(^) = limE^[(P^x,P^B(^x))]. 
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Since x = P' K x + F h K x, 

(T l K x,T l K B(x,x)) = (F l K x,B(x,F l K x)) + (F> K x,B(x,F h K x)} 

= (T l K x,B(x,P h K x)) 

as sum of a finite number of terms (so we can use the 

anti-symmetric property of the non-linear term without convergence issues). 
In conclusion, 
(4.9) 

l( J L^) = limE^[(Pj i: x, J B(jc,P|x)}] = limE^[ £ (x m -xj)(m-X[ 



l+m=k 

I 

\mU>K 



Following Frisch [15, Section 6.2], the last term we have obtained in the 
formula above is the energy flux through wave-number K and represents 
then energy transferred form the scales up to K to smaller scales. 

From the previous section we know that i (/!*) > 0, this is a consequence 
of property [M3] of Definition 2.4. From a mathematical point of view, ex- 
istence of invariant measures with I (/I*) > would be an evidence for loss 
of regularity and, in turn, for blow-up. From a physical point of view, the 
energy flux through wave-numbers should converge to zero - hence, again 
we would expect l(/l*) = - as the energy should flow through modes es- 
sentially only in the inertial range (we refer again to Frisch [15]). 

Proposition 4.11. We have 

1. the map /l* i— > e(jU*) has a smallest element in ( solution of largest 
mean inertial dissipation), 

2. if 



lim sup ^[{\x\ 2 v -R)l {K>R} ]=0 



then there is /!* G J? such that £(/!*) = inf^ e y m £(/!*) (solution of 
smallest mean inertial dissipation). 

Proof. The first part follows easily as J 1 is compact (see Remark 4.2) and 
\x — > e(jii) is lower semi-continuous for the topology with respect to which 
J* is compact. As it regards the second part, we know by Corollary 3.2 that 
J'm is relatively compact on Q,(V). Hence, if M is the largest value attained 
by £ on J? m and (/l, 3 )„ e N is a maximising sequence, say e(fl n ) > M — ^, by 
compactness there is /i* such that, up to a sub-sequence, /i„ — > /l^. Now 

e R (ji n )=W^[\x\ll m 2 m ]+RfJ fl [\x\l>R] 

>M-~- sup E»[(\x\ 2 v -R)l { , x a >R} ], 

where £r(/i) = E^[\x\yAR]. As n | °°, by continuity £r(jU„) — ► £r(p^), so 
£(jUoo) >£^(/ioo) >M- sup / Wvl/u2 >i? } /x(<±c). 
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As R — > oo, it follows that > M, hence £(jUoo) = M. □ 

We have not been able yet to prove the condition given in item (2) of 
previous proposition. 

We also remark that such measures of largest and smallest mean iner- 
tial dissipation may not be unique, as both functionals £(•) and t(-) are 
translations-invariant (see Section 4.2). 



Appendix A. Analysis of the mildly regular approximated 

PROBLEM 

Let R > 1 and let %r : [0,°°] — > [0, 1] be a non-increasing C°° function 
such that %r = \ on [0, jR], %r = on [2/?,°°) and there is c > such that 
\Xr\ < f (see the picture). 



\R 2R 
Figure 1 . The cut-off function %r 



Given a value £0 G (0, |], we consider the following problem in ^ eo , 

(A.l) \ £ o 
divJ?*'*) =0. 



Let T (e °' s) : Q — »• [0, 00) be defined as 



T^ R \(o)=mf{t>0:\co(t)\l EQ > 3 -R} 

(and T (£ ° ■ R \co) = 00 if the above set is empty). The main aim of this section 
is to analyse the solutions to the above problems and their connections to 
the original Navier-Stokes equations (1.2). 

Before turning to the results on the regularised problem (A. 1), we remark 
that in the proof of all results of this section we shall use the splitting u (e "- R) = 
y(Eo,R) _|_ ^ w here z solves the following linear Stokes problem 



(A.2) 



dz+Azdt = 0,2 dW, 



k z(0) = 0, 

and so v^ - R) solves the following equation with random coefficients 
(A.3) ^ £b - r) + vA^- r) +Xr(\^ £o - r) \1 )B(rf £ °- R \rf E ° R) )=0. 
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A.l. The weak-strong uniqueness principle. We first extend the weak- 
strong uniqueness principle stated in [12, Theorem 5.4] to the above prob- 
lem (A.l). This is the content of the following result. 

Theorem A.l. Assume condition [A2] of Assumption 1.1 and let £o G (0, ^] 
with £q < 2(Xq. Then, for every x G n/ £o equation (A.l) has a unique martin- 
gale solution Px , with 

(A.4) Pi^[C([0,=o);^)] = l. 

Moreover, the following statements hold. 

1. (weak-strong uniqueness) On the interval [0, t {£ "- R} ], the probability 
measure °' coincides with any martingale solution P x of the orig- 
inal stochastic Navier-Stokes equations (1.2), namely 

E P ^ R) [fl>(6)l {T (*,*)> f }] = E ft [<pfe)l { *«,*)> t} ] 

for every t > and every bounded measurable <p : H — >■ R. 

2. (blow-up estimate) There are cq > 0, c\ > and C2 > 0, depending 
only on £q, such that for every x G with |x| 2 £ < j, 

(A.5) pfa*)[ x M <§]< Cie - C 2x 

i_ 

for every 5 < cqR 2£ o . 

Proof. The proof is developed in four steps, which are contained in the fol- 
lowing subsections. More precisely, in the first step we prove existence of 
solutions for problem (A.l), while in the second step we prove uniqueness. 
The weak-strong uniqueness principle is then proved in the third step and 
the blow-up estimate (A.5) is given in the fourth step. 

Step 1 : Existence. We only show the key estimate for existence. Let z be 
the solution to the linear Stokes problem (A. 2) and consider v^ ,R) as above. 
The usual energy estimate provides (here we use u = w (£o,R) and v = v^*) 
for brevity) 

+2vl|vll|. ab < 2^(|«q^|)<A*-iB(i*;ii),A5+*v>a 
<cXR(M\\)H^\A e ^u\ 2 
(A.6) <^(l^^l)IMI^(|A e( ^z| 2 +N^|M|^ 2eo ) 

<v||v1|^ + c(|A ^| 4 +|4^+^), 

where we have used interpolation inequalities and Lemma 2.2, with a = Eq. 
Notice that, by the choice of Co with respect to Oq, |A e ^z|^ has exponen- 
tial moments. 

In order to show (A.4), we show an a-priori estimate for the derivative 
in time Ea ' R) in L 2 (0, T;D(A-^- £ ^)), for all T > 0. The continuity of 
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jjfa>,«) then follows from this fact and continuity of z- The a-priori estimate 
follows by multiplying the equations by ^A 2£6_ 2v (£|,! ' i) , 

where we have used the same estimates as in (A. 6) (and again u = u i£o ' RS > and 
v = y^*) for brevity). 

Step 2: Uniqueness. Let u\, u,2 be two solutions of (A.l) starting at the 
same initial condition and set w = u\ — u 2 - The new process w solves the 
following equation with random coefficients, 

w + vAw = Xr(\u2\% )B(u 2 ,u 2 ) - Xr{\u\\\, f )B(u h ui) 

= [Xr(\u 2 \1 £() ) -Xr(\ui\1 £o )]B(u u u 2 ) + Xr(\u 2 \1 £o )B(w,u 2 ) 



+ Xr(\ui\1)B(ui,w) 



and so 

d 



(lt Ml + ^\Ml = ^[XR(\u2\l £ J-XR(\ui\l £Q )](w,B(u u u 2 )) H 
+ 2xr{\u 2 \IJ(w,B{w,u 2 )) h 



=0+0 



Next, we estimate the two terms. In order to estimate the first term, we first 
remark that 

\Xr(\u2\1 £q ) ~ XR(\m\l £Q )\ < -^M^[l [ o i 2i ? ](|«l£ £() ) + l[0,2i?](l"2|^)]. 



By using the above inequality, Lemma 2.1 (with OL\ = j — Co, Ct 2 = and 
0C3 = | + £o) and interpolation and Young's inequalities, we get 5 

[3< \XR(\u2\l^)-XR(\ui\\)\\(B{u h w),U2)H\ 

< ^\^\w\^ \w\ v (\A^ £o ui\ H + \A^ £ °u 2 \ H ) 



^o 1 

1 -2e 0| ,|+2cb/| A l- 



<c\w\l HHy (i A2 ~^"i|H + l A3 ~^"2|ff) 



V 



-2 i 

The inequality has to be slightly modified if 6q = j - In such a case we use Lemma 2. 1 



with ai = j, a.2 = 0, a3 = ^, 
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For the second term, we use Lemma 2.1 with a\ = i + Co, 0C2 = and 
(Z3 = 5 — Co, and interpolation and Young's inequalities, 

\2\ = 2Xr{\u2\\)\{u 2 ,B{w,w)) h 

< 2cXr{\u 2 \1 Eq )\u 2 \v £o \w\v \A^ £ °w\ h 

(A.7) < 2cR\w\u Q \ W \ 2 v {l - £i)) 

1, ,? J -, ,? 

< -\w\y + cR £ o \w\fj. 

Finally, Gronwall's lemma implies that w = 0, since w(0) = 0. 

Step 3: WeaA-sfrong uniqueness The proof works exactly as in [12, Theo- 
rem 5.12], we give a short account for the sake of completeness. The proof 
is developed in the following steps. 

1. The energy balance of w = m (£ ° r) — ^ , given by 

g t = \w{t)\ 2 H + 2v f \w(s)$ds+ 
Jo 

+ 2 f^XR^X^B^ 

is an a. s. super-martingale under P x . 

2. t (Eo S) is a stopping time with respect to the filtration {^§ t )t>Q- 

3. The stopped process (^ AT (£„^))f>o is again an a. s. super-martingale. 

4. The previous step implies the conclusion. 

All the above steps can be carried out exactly as in the proof of Theorem 
5. 12 of [12] (the key point is that ifo - R) is continuous in time with values in 
^e with probability one). The only difference is in the last step, where the 
estimate of the non-linearity needs to be replaced by the following estimate, 

(w,B(^ R \^) -B(Z,S)) H = -(^> R \B(w,w)) H . 

Finally, the above estimate can be obtained as in (A.7). 

Step 4: The blow-up estimate Fix x G ^ eo with \x\ 2 v < § and 8 > 0. Set 

&§ = sup se [ 31 |A e ( e °)z| 2 , then, by slightly modifying inequality (A. 6), we 
get 

j t \^ R) \\ < c(\A e ^z\ 4 H + \^ R) C £ ;) < c( 1 + 2 5 + |^%)^, 

where v< £o s) = 5**° - R) — z has been defined in (A. 3). Hence, if we set (p(t) = 
^ + 0g + l^^'^'l^ , we end up with a differential inequality that, once 
solved, gives 

<p(t) < <p(0)(l -c8(p(0)^y 2£Q 
From this, it is easy to show that there is a suitable constant cq = cq(£q) > 
such that |S tEo - R) (s) | 2 < jR for every s < 8, when &§ < § and 8 < cqR^ 1 ^ . 
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In conclusion, if 8 < cqR~^ and \x\ 2 , < f , then 

5 <| 1^(01^ <^ for?<5 =^ t (£ »' r) >5 

and so 

pfo>>*)[ T C*.«) < §] < [® g > *] < cie" C2 ^ 

with constants ci = ci(£o) > and C2 = C2(£o) > depending only on £o 
and the last inequality follows easily as in Proposition 15 of [ ] (which in 
turns is a consequence of Proposition 2.16 of Da Prato & Zabczyk [ ]). □ 

A. 2. Moments of norms of stronger regularity. The proof of the follow- 
ing theorem is based on an inequality taken from Temam [22, Section 4.3, 
Part I] (see also Odasso [18]). 

Theorem A.2. Under condition [A2] of Assumption 1.1, for every £q 6 (0, |], 
with Eq < 2(Xo, there are 5 = 5(£q, OCq) > and y = y(£o, OCq) > such that 



(A.8) E 



7V^ R) | 2 i ds\ <C[\+t + \x\ 2 H ], 

1/0 a J 



where w (Eo,s) is the solution to problem (A. 1 ) starting atxE and the value 
ofC is independent of both x and R. 

Proof If OCo < \, we choose £o < 2«o (such condition is useless for all 
other values of OCq). The noise is not regular enough to let us work directly 
on w (£o R) , so we rely, as in the proof of the previous theorem, on yte>' fi ) — 



g(«b.*) _ Zm L et p — 1 _ £o (the value of p could be slightly improved, but it 
is beyond our needs) and compute 



"rn^n2 \-m ~ m \ „ n ~a ^ {A2 +2eo v,B(u,iP)} H 
-[(1 + 1^) P \ = 2P 2 -IpXRil^) : , t , 

where we have set w = ^ £ ° fi) and v = v^^'. The non-linear term can be 
estimated as in (A. 6) to get 

(A«*+!v,A*-ifl(5f,2)> fl < ^l|v|| 2 , £o +c(|A e ( e «)z|4 +|v| 2 ^) 

and so 

^\ + . ds<l+c f\A 6 ^z\ 4 d S + c f\v\l ds. 
Jo (l + |vl^ o ) p+1 Jo Jo £() 

The term in z is plain (see for example Da Prato & Zabzcyk [ ]), while the 
term in \v\^ can be estimated by means of the energy inequality [M3] of 
Definition 2.4. Finally, in order to prove (A.8), we use again the energy 
balance, since by Young's inequality 
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if 7 is chosen properly, depending on p. 

Next, if GCo > k' fi x a e (?) ^b) an d £o G (0, 5] and choose p > (whose 
value will be fixed in dependence of a). We apply Ito formula to the func- 
tion (1 + \A a ^- R \t)\ 2 H )-P, to get 

1 1 

(1 + \A a u(t)\ 2 H )P ~ (1 + |A«a:|^)p ~ 

io (l + |A« M |^)/'+ 1 7o (l + \A a u\ 2 H )P +1 



and we have set again u = m (£o,s) . The non-linear part is estimated with 
Lemma 2.2, interpolation and Young's inequalities, 

(A a+ ?u,A a -W(u,U)) H < C\A a+ ?u\ H \A e{a -^u\ 2 H 

< C\A a +-2u\ H \A a ^l a+ ?\A a+l 2u\l- 2a 

1 1 „ o 4a+l 

< ^\A a+ 2u\ 2 H+ c\A a u\^~ l . 

If « < \ + Co, a + j > 0(oco) an d one already knows that some power of 
|A a w| has finite moment, then one can proceed as in the previous case (Xq < 
\. Otherwise, as in Temam [2 ], one can iterate the same procedure using 
a — 5 instead of a, until the above conditions are satisfied. □ 

A. 3. An estimate of the return time to a ball. The aim of this section 
is to verify that the probability of hitting a ball (in a smooth norm) can be 
uniformly bounded from below for all initial condition in a given ball. 

Lemma A.3. Assume condition [A3] from Assumption 1.1. Then one can 
choose Co G (0, \] with £$ < 2(Xq such that there are Q' < 0Cq + \, 9" > 0((Xq) 
and a suitable constant c > 0, and the following statement holds. 

For every R > 1 there are values Tq = Tq(R) and K = K(R) such that for 
every x G <V £q , 



x\l <R, 

ifl',C.M2 ^ n b. r_, J 1 ^ J 0) 



w su PreM |A« z (,)\i,<R =► m{ |A ^ sw , (Ii)| 2 <jri 

[T <cR 

where z is the solution to the linear problem (A. 2). 

Proof. We choose £0 = 5 and we set, for brevity, u = w (£o ' 3 ^ and v = u — z. 
The first part of statement [B] follows as in the proof of (A. 5), if the constant 
c is chosen accordingly. So, for every t G [0, 7b], we know that \u(t) \ v < 3R. 
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In particular, using (A. 6) and the second statement of [A], it follows that 
there is a constant Kq = Kq(R) such that 

(A.9) sup \v(t)\ 2 v + [ T °\Av(s)\ 2 H ds<Ko(R). 

[0,2b] J ° 

We next prove the second statement of [B]. 

Step 1. We first consider Oq G (5,3] and we choose 5 G (0, |) so that 
^ < 8 < 8 + \ (this condition ensures that 0(ctb) < 3 + 5 and 0(5 + \) < 
(Xq + 5 ). In the case a = \ we simply choose a value 5 G Q, |). 

Step 2. For all CO £ £1 satisfying [A], there is to = to(co) G (0,Tq) such 
that \Av(to)jj < 2K . Indeed, from (A.9) it follows that the set {t e (0, T ) : 
\Av{t) 2 H < 2Kq} has Lebesgue measure at least -|, and in particular is not 
empty. 

Step 3. There is K\ = K\{R) such that for all co E Q. satisfying [A], 
\A? +s u(To) fa < K\. We use Lemma 2.2 (with a = 8 + |), interpolation 
of D(A e ( s+ ^) between V and D(A l+s ) and Young's inequality to obtain 
the following estimate, 

||aH^ + 2v|A 1 +^ = 2(a i + 2 5 ~ jB( ~^ )) 

< c\A l+5 v\ H \A e ^u\ 2 H 

<v\A l+s y H + C(\A e ^z\ 4 H + M 6 v +sS ). 

Since v is bounded in V, the claim easily follows. In the special case Cto = j 
one can proceed analogously. 

Step 4. We choose then 9' = 9(8 + \) and 9" = 8 + \ and the second 
statement of [B] follows. 

Step 5. If Ofo > i we iterate the above procedure as in the proof of Lemma 
C.l of [12], using the two inequalities 

[il] if m E N, m > 1, there are an integer p m and C m > such that 

j t \A*v\ 2 H + v\A^ 2 H < C m {\ + \v\ v + \A%z\ H ) Pm , 
[i2] if K > I and /3 G [0, |), there are C k r > and ap > such that 

||A^vl| + v|A^ + M|<C M [|A^z|l + (|A' f+ ^||A' c v1^ 

+ \A^z\ 2 H )\A^v\% 
whose proof can be found in the proof of that same lemma. □ 

References 

[1] D. Blomker, F. Flandoli, and M. Romito, Markovianity and er- 
godicity for a surface growth PDE (2006), preprint, available at 
http : //arXiv . org/abs/math . PR/061 1021. 

[2] G. Da Prato and A. Debussche, Ergodicity for the 3D stochastic Navier-Stokes equa- 
tions, J. Math. Pures Appl. (9) 82 (2003), no. 8, 877-947. 



EQUILIBRIUM STATES OF MARKOV SOLUTIONS TO NAVIER-STOKES 



31 



[3] G. Da Prato and J. Zabczyk, Stochastic equations in infinite dimensions, Encyclope- 
dia of Mathematics and its Applications, vol. 44, Cambridge University Press, Cam- 
bridge, 1992. 

[4] G. Da Prato and J. Zabczyk, Ergodicity for infinite-dimensional systems, London 
Mathematical Society Lecture Note Series, vol. 229, Cambridge University Press, 
Cambridge, 1996. 

[5] A. Debussche and C. Odasso, Markov solutions for the 3D stochastic Navier-Stokes 
equations with state dependent noise, J. Evol. Equ. 6 (2006), no. 2, 305-324. 

[6] C. Dellacherie and P. -A. Meyer, Probabilities and potential. B, North-Holland Math- 
ematics Studies, vol. 72, North-Holland Publishing Co., Amsterdam, 1982. Theory 
of martingales; Translated from the French by J. P. Wilson. 

[7] J. Duchon and R. Robert, Inertial energy dissipation for weak solutions of incom- 
pressible Euler and Navier-Stokes equations, Nonlinearity 13 (2000), no. 1, 249-255. 

[8] F. Flandoli, On the method of Da Prato and Debussche for the 3D stochastic Navier 
Stokes equations, J. Evol. Equ. 6 (2006), no. 2, 269-286. 

[9] F. Flandoli and D. Gatarek, Martingale and stationary solutions for stochastic 
Navier-Stokes equations, Probab. Theory Related Fields 102 (1995), no. 3, 367-391. 

[10] F. Flandoli, M. Gubinelli, M. Hairer, and M. Romito, Rigorous re- 
marks about scaling laws in turbulent fluids (2005), preprint, available at 
http : //arXiv . org/abs/math-ph/0507044. 

[11] F. Flandoli and M. Romito, Markov selections and their regularity for the three- 
dimensional stochastic Navier-Stokes equations, C. R. Math. Acad. Sci. Paris 343 
(2006), no. 1,47-50. 

[12] F. Flandoli and M. Romito, Markov selections for the three-dimensional stochastic 
Navier-Stokes equations (2007), to appear on Probab. Theory Related Fields. 

[13] F. Flandoli and M. Romito, Regularity of transition semigroups associated to a 3D 
stochastic Navier-Stokes equation, Stochastic Differential Equations: Theory and Ap- 
plication (P. Baxendale H. and S. Lototski V., eds.), Interdisciplinary Mathematical 
Sciences, vol. 2, World Scientific, Singapore, 2007. 

[14] F. Flandoli and M. Romito, 2007, in preparation. 

[15] U. Frisch, Turbulence: The legacy of A. N. Kolmogorov, Cambridge University Press, 
Cambridge, 1995. 

[16] B. Goldys and B. Maslowski, Exponential ergodicity for stochastic Burgers and 2D 
Navier-Stokes equations, J. Funct. Anal. 226 (2005), no. 1, 230-255. 

[17] S. P. Meyn and R. L. Tweedie, Markov chains and stochastic stability, Communica- 
tions and Control Engineering Series, Springer- Verlag London Ltd., London, 1993. 

[18] C. Odasso, Spatial smoothness of the stationary solutions of the 3D Navier-Stokes 
equations, Electron. J. Probab. 11 (2006), no. 27, 686-699 (electronic). 

[19] C. Odasso, Exponential mixing for the 3D stochastic Navier-Stokes equations, 
Comm. Math. Phys. 270 (2007), no. 1, 109-139. 

[20] M. Romito, Existence of martingale and stationary suitable weak solu- 
tions for a stochastic Navier-Stokes system (2001), preprint, available at 
http : //arXiv . org/abs/math . PR/0609318. 

[21] D. W. Stroock and S. R. S. Varadhan, Multidimensional diffusion processes, 
Grundlehren der Mathematischen Wissenschaften [Fundamental Principles of Math- 
ematical Sciences], vol. 233, Springer- Verlag, Berlin, 1979. 

[22] R. Temam, Navier-Stokes equations and nonlinear functional analysis, 2nd ed., 
CBMS-NSF Regional Conference Series in Applied Mathematics, vol. 66, Society 
for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, 1995. 



32 



M. ROMITO 



dlpartimento di matematica, universita di flrenze, vlale morgagni 
611 a, 50134 Firenze, Italia 

E-mail address: romitoOmath . unif i . it 



